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ABSTRACT 



We use mid-IR images from the Spitzer Space Telescope Cygnus X Legacy 
Survey to search for stellar bowshocks, a signature of early type "runaway" stars 
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with high space velocities. We identify ten arc-shaped nebulae containing cen- 
trally located stars as candidate bowshocks. New spectroscopic observations of 
five stars show that all are late O to early B dwarfs, while one is a previously 
classified BO. 2 giant. These stars have moderate radial velocities, differing by 
AV <10 km s _1 from members of the Cygnus OB2 Association. The spectral 
energy distributions of the other four stars are consistent with late O to early B 
dwarfs at the nominal ~1.6 kpc distance of Cyg OB2. Our morphologically se- 
lected sample of bowshock candidates encompasses diverse physical phenomena. 
Three of the stars appear to be pre-main-sequence objects on the basis of rising 
SEDs in the mid-IR, and their nebulae may be photon-dominated regions (PDRs) 
illuminated by the central star but shaped by external sources such as winds from 
Cyg OB2. Four objects have ambiguous classification. These may be partial dust 
shells or bubbles. We conclude that three of the objects are probable bowshocks, 
based on their morphological similarity to analytic prescriptions. Their nebular 
morphologies reveal no systematic pattern of orientations that might indicate ei- 
ther a population of stars ejected from or large-scale hydrodynamic outflows from 
Cyg OB2. The fraction of runaways among OB stars near Cyg OB2 identified 
either by radial velocity or bowshock techniques is ~0.5%, much smaller than the 
~8% estimated among field OB stars. We discuss possible reasons for this dif- 
ference. We also obtained a heliocentric radial velocity for the previously known 
bowshock star, BD+43°3654, of -66.2 ± 9.4 km s _1 , solidifying its runaway sta- 
tus and implying a space velocity of 77±10 km s _1 . We use the principles of 
momentum-driven bowshocks in conjunction with the observed sizes, bowshock 
luminosities & spectral energy distributions, and dust/PAH emission models to 
arrive at a novel method for estimating stellar mass loss rates. Derived mass loss 
rates range between 1CT 7 and fewxlCT 6 M yr _1 for the three 05V - ~B2V 
stars identified as generating bowshocks. These values are at the upper range 
of, but broadly consistent with, estimates from other methods. We calculate a 
relatively large mass loss rate of 160xl0~ 6 M yr" 1 for 04If star BD+43°3654 
using the same method. 

Subject headings: stars: early-type - stars: individual: BD+43°3654, HD 195229, 
GSC 03161-01188-Galaxy: open clusters and associations: individual: Cygnus 
0B2 
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Introduction 



Most OB stars lie within the associations in which they were born. However, a significant 
fraction (17%-50%) having high (30-200 k m s" 1 ) space velocities wander far from such 
associations and are known as " runaways" (IBlaauwl Il96ll ; IStond Il982l ; IConti et al.l 119771 ; 
Gies &: Bolton! Il986l ). Often, these stars' proper motions allow them to be traced back to a 
known OB association in which they presumably formed. The speed necessary for a star to 
be classified as a runawa y varies from 30-40 km s _1 depending on author. We adopt the 30 
km s -1 criterion used by lGies fc Boltonl (119861 ). 



Blaauwl (119611 ) suggested that runaway stars are ejected from OB associations when one 



member of a binary system sheds a large fraction of its mass, as in a type II supernova. In 
this case (known as the binary-supernova scenario), the survivin g comp anion is released with 
a velocity comparable to its previous orbital velocity. IBlaauwl ( 19611) cited the low binary 
frequency of runaway stars as evidence for this model, though I Gies fc Boltonl (119861 ) found 
a few runaway binaries in their survey. Whether and how frequently supernovae in binary 
systems produce an unbound, high-velocity O star is an active line of investigation. 



Gies fc Boltonl (119861 ) proposed that runaway stars are ejected during close encounters 



between binary systems in OB associations. This model, known as the dynamical-ejection 
scenario, explains the existence of binary runaway systems, since hard encounters among 
three or more stars can produc e both single runaway stars and binary runaway systems. 
Leonard fc Duncan Jl988l . Il990l ) performed numerical simulations of young star clusters, 
showing that dynamical encounters can readily account for the number of observed runaways 
and that runaway binary systems are possible. 



Hoogerwerf. de Bruijne. fc de Zeeuwl (120001 ) used precision proper-motion measurements 



of nearby runaways to find evidence for both binary-supernova and dynamical ejection mech- 
anisms. They concluded that the runaway ( Oph likely resulted from the dissociation of a 
binary system wherein the other member became the pulsar PSR J1932+1059. They also 
concluded that the runaways AE Aurigae and /i Columbae and the binary system i Orionis 
were produced in a binary-binary encounter. 

Traditionally, runaway star s reveal themselves t hrough their high proper motions o r 
unusually large radial velocities (IGies &: Boltonlll986l ; iMdzinarishvili fc Chargeishvilil 120051 ) . 
Identification using radial velocities requires high-precision spectra, and this technique is 
sensitive only to runaways with large projected radial velocity components. Identification 
on the basis of hi gh proper motions req uires precise astrometric data, such as those made 
by the Hipparcos (jPerryman et al.lll997l ) mission, along with secure distance measurements. 
This technique is only sensitive to runaways with large tangential motions. 
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Runaway O stars have also been identified by detecting the tell-tale "bowshocks" pro- 
duced when the stellar winds of stars traveling supersonically impinge upon the surround- 
ing interstellar medium. Bowshocks appear as symmetric arc-shaped features with apsides 
oriented in the direction of the stars' motions and at distances from the star determined 



by momentum balance betwee n th e wind and the ambient med ium. Ivan Buren fc McCray 



( 119881 ) . Ivan Buren et al.l (119951 ). and lNoriega-Crespo et al.l (119971 ) used Infrared Astronomical 



Satellite (IRAS) images to locate arc-shaped features associated with high-velocity O stars. 
These studies noted that a variety of phenomena can produce bowshock-like morphologies, 
including partial stellar wind bubbles and dust shells, H II regions with density gradients, 
and genuine high-velocity stars. Indeed, we conclude that our sample is comprised of a mix- 
ture of such objects. Several more recent studies have used Ha surveys, Midcourse Space 



and associate them with high-velocity stars ( 


Brown & Bomans 


2005; 


Comeron & Pasauali 


2007; 


Gvaramadze & Bomans 


2008; 


Povich et al. 


2008). Bowshocks are also seen in envi- 



ronments such as the Orion Nebula where winds and j ets from young en counter the O-star 
winds that power the nebula (e.g., the HST images of iBally et all (boooh ). 



In this paper we use infrared images from the SST Cygnus-X Legacy Survey (IHora et al. 



20071 ) to identify bowshock candidates in the vicinity of one of the Galaxy's richest OB 
associations. The association Cygnus OB2, at t he heart of the Cygnus X r e gion, contains 



over 100 OB stars, including an Q3If and Q4If (IMassey &; Thompson! 1 1991 



Hanson 2003 



Comeron et al.l 12002; IComeron fc Pasaua 



20071 ). The spectroscopic OB survey of 142 stars 



in Cygnus OB2 by iKiminki et al.l (120071 ) did not detect any runaway stars on the basis 
of radial velocities. However, in a region as rich in young massive stars as Cygnus OB2, 
one would also expe c t to find stars ejected with large motions in the plane of the sky. 
Comeron fc Pasqualil (120071 ) implicate Cygnus OB2 as the origin of one such star, the 04If 
runaway BD+43°3654, using the morphology of the bowshock seen in MSX images. Given 
the superior sensitivity and angular resolution of SST, a more complete search may now 
be conducted for runaway stars hosting bowshocks in this region. Herein, we describe the 
discovery of ten candidate bowshocks near Cygnus OB2, and we report followup optical 
spectroscopy that allows us to determine the spectral types and radial velocities for a subset 
of the central stars that power each candidate. We refer to the bowshock candidates and 
their central stars separately as Bowshocks (BS) 1-10 and Stars 1-10. We also investigate 
a novel method of determining the mass loss rates of these stars using the principles of 
momentum-driven bowshocks. 



Modern distance estim ates to Cy g OB2 range f rom 1.7 kpc (IMassey fc Thompsonlll991 



Torres-Dodgen et al. 


1991 


) to 1.4 kpc ( 


Hanson 


2003) 



1 pc corresponds to 130" on the sky. 
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Data 



2.1. Spitzer Space Telescope 



We retrieved data from the SST archive obtained as part of the Spitzer Legacy Survey 
of the Cygnus-X Complex (jHora et al.l 120071) using data fro m the Infrared Array Camera 
(IRAC) at 3.6, 4.5, 5.8, and 8.0 microns ([Fazio et all 12004 ) an d the Multi-band Infrared 
Photometer for Spitzer (MIPS) bandpasses at 24 &: 70 microns (IRieke et al.l 120041 1 . One of 
us (HAK) visually examined mosaiced single-band and multi-color images of the ~24 square 
degree survey area to identify bowshock candidates on the basis of symmetric arc-shaped 
structures less than a few arcminutes in size enclosing a symmetrically placed stellar source. 
The complex ISM structure in this region, coupled with the large dynamic range of features, 
renders bowshock identification a subjective and imprecise task. In practice, the 8-micron 
band (sensitive to emission from PAH features excited by non-ionizing photons) and the 
24-micron band (sensitive to emission from warm dust) were the principal images used to 
locate candidate bowshocks. While a plethora of arc-like structures appear in both bands, 
very few exhibited a high degree of symmetry and included a point source near the apsis. 
In total, ten objects were deemed to be bowshock candidates and were retained for further 
investigation. Table UJ lists the equatorial &; Galactic coordinates and 2MASS JHK & SST 
IRAC photometry of the central stars for each object. The final column provides common 
stellar cross identifications from other works, where applicable. 

Figured] is a 3-color image ~ 2° square centered on Cygnus OB2 (£ = 80.3°, b = +1.0°) 
showing the location of the first seven candidate bowshocks. Blue, green, and red represent 
the 4.5, 8.0, and 24 micron images, respectively. Objects 8-10 have Galactic latitudes b < 
0.15° and lie off the lower edge of Figure [TJ Figures [2] through [TTJ show enlarged views of 
each bowshock candidate using the same color scheme as Figure [TJ We discuss each object 
in detail in §3. 



2.2. WIRO 

We obtained optical spectra for several of the stellar sources located near the apsides of 
candidate bowshocks. We observed five of the stars (Stars 1, 2, 3, 5, and 8) at the Wyoming 
Infrared Observatory's (WIRO) 2.3 m telescope equipped with the Longslit spectrograph on 
2008 June 23-30, 2009 May 20, and 2009 September 16. The 1800 1 mm" 1 grating was used 
in first order. The spectral coverage was 5210-6680 A. Exposure time varied from 360 s for 
Star 3 (HD 195229) to 3600 s for Stars 2 and 5, yielding signal-to-noise ratios ranging from 
50:1 for Star 5 to 450:1 for Star 2. Copper-Argon lamp exposures were taken after each 
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stellar spectrum for wavelength calibration. 

The spectra were reduced, extracted, wavelength calibrated, and Doppler corrected to 
the heliocentric frame using; standard IRAF0 tasks in the KPNOSLIT package. The absolute 
radial velocities were checked against the standards HD161096, HD 161797, HD131156, 
HD171391, and HD188512. Table[T]lists the heliocentric radial velocities, which are typically 
accurate to 6 km s" 1 rms. 



2.3. WIYN 

We also obtained spectra of BD+43°3654 and HD195229 during an observing run at 
the WIYN@ 3.5 m telescope on 2008 June 11-16 using the Hydra multispectrograph. We 
used the Red camera and 3" blue fiber cables with the 1200 1 mm -1 grating in second 
order. Exposure times were 600 s and the spectral coverage was 3820-4500 A. Copper- 
Argon lamp exposures were taken to wavelength calibrate the spectra. The spectra were 
reduced, extracted, wavelength calibrated, and Doppler corrected to the heliocentric frame 
using standard IRAF tasks in the TWODSPEC and RVCOR packages. Radial velocities 
of standards HD 131156, HD146233, HD161096, HD161797, and HD 171391 from this run 
agreed with published values to ±2 km s _1 . 



3. Analysis 
3.1. Bowshock Identification and Images 



Figure [T] shows a portion of the Cygnus-X region surrounding its most massive and 
energetic constituent, Cygnus OB2. This region encompasses a complex array of neb- 
ulae and star- forming regions located at various distances along the line of sight, with 
most of the activity located in the local or "Orion" spiral arm at a distance of 1-2 kpc. 
Odenwald fc Schwartzl Jl993h pr esent a schematic 3-dimensional representation of the Cygnus- 
X region, while ISchneider et al.l (120061 ) provide velocity-resolved molecular maps of Cygnus-X 
and individual star-forming objects nearby. One apparent feature of this region is the rel- 



1 IRAF is distributed by the National Optical Astronomy Observatories, which are operated by the As- 
sociation of Universities for Research in Astronomy, Inc., under cooperative agreement with the National 
Science Foundation. 

2 The WIYN Observatory is a joint facility of the University of Wisconsin-Madison, Indiana University, 
Yale University, and the National Optical Astronomy Observatories. 
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ative absence of dust (red) and PAH emission (green) at the location of Cyg OB2 itself, 
except at the southern edge. This apparent "hole" is consistent with a cavity cleared by 
the winds of >100 evolving massive stars. Some of the prominent star-forming regions la- 
beled in Figure [T] clearly constitute the heads of gaseous pillars, most of which point toward 
Cyg OB2, the dominant source of lumi nous and mechanical energy in this vicinit y. These 
pillars are reminiscen t of those in M 16 Jlndebetouw et al JbooH Ikester et alJll996h . RCS 49 
( jWhitney et al.l 120041 ) and other massive star-forming regions. 

The seven bowshock candidates labeled in Figure Q] lie around the periphery of Cyg 
OB2, none falling within the main concentration of Cyg OB2 stars but all within ~5-8 
pc of the canonical Cyg OB2 center. The orientations of the features, seen more clearly 
in Figures [2] through [TTJ appear random, having no preferential alignment with respect to 
Cyg OB2. Generically, bowshock orientation reveals the relative motion between the star and 
surrounding matter. We might have expected, if the stars were all runaways from Cyg OB2, 
that the bowshocks would point away from the Association in the direction of the stars' 
motions. On the other hand, if the hydrodynamics of material surrounding Cyg OB2 were 
dominated by a hot outflow from an over-pressured region, we might have expected many of 
the bowshocks to point back toward the source of mechanic al energy, as do three bowshocks 
near M 17 and at least one near RCW 49 (IPovich et al.ll2008l ). The absence of either signature 
suggests that 1) the central stars driving the bowshocks, if they are runaways, have different 
origins, and 2) there is no large-scale supersonic wind emanating from Cyg OB2. Only 
candidates 2 and 4 have orientations consistent with their stars being runaways having an 
origin within Cyg OB2. 



3.2. Nebular and Stellar SEDs 

The colors of the candidate bowshock nebulae , best seen in the individual zoomed 
images in Figures [2] through [TTJ at first glance reveal two classes of objects. Candidates 
1,3,4,7 are seen exclusively at 24 microns and longer wavelengths while the rest exhibit strong 
emission in the IRAC bandpasses with some additional contribution at longer wavelengths. 

We performed aperture photometry of each bowshock candidate in each of the six 
bandpasses using the SSC post-basic calibrated data mosaic images by assigning irregular 
crescent-shaped polygonal apertures defined visually on the [8.0] and/or [24] images. We also 
manually defined background regions surrounding each bowshock. Because these objects lie 
in complex regions of diffuse emission, care was taken to select multiple background regions 
judged to be representative of the local background away from the nebula but within several 
arcminutes. We applied aperture corrections, as recommended in the IRAC and MIPS cali- 
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bration web pageg§, using values appropriate to a circular aperture of equivalent area to the 
polygonal apertures. These corrections are less than 1.3 in all cases and are sometimes less 
than unity for the IRAC bands. Given the imprecision of the aperture corrections and the 
high and variable background levels, we adopt a minimum uncertainty of 20% in all bands. 
Table [2] gives background-subtracted aperture fluxes in mJy for all of the nebulae. In some 
cases, only la upper limits are given, mostly in the IRAC bands where some objects are not 
detected. 

For each bowshock candidate we were able to identify a symmetrically placed point 
source, presumed to be the energizing star, behind the apsis. Figure [12] is a 2MASS JHK 
color-color diagram showing each of the ten stars (diamonds) associated with the bowshocks. 
The solid and dashed curves (asterisks and crosses, respectively) show a fiducial main se- 
quence and supergiant sequence. A solid line illustrates the reddening vector for Ay = 5 
mag. The majority of stars lie in a close group near H-K=0.3, J-H=0.6, consistent with 
early-type stars seen behin d 4-7 magnitudes of yisual extinction. This range is similar to 



other OB stars in Cyg OB2 (IMassey fc Thompson! 1 19911 ). suggesting that they lie at a similar 



distance. Stars 4 and 6 are much redder, suggesting 12-15 magnitudes of visual extinction, 
under the assumption that these are also early-type stars. The K-band magnitudes for these 
two stars are 2-3 magnitudes fainter than the other targets, indicating either a larger dis- 
tance or the effect of localized regions of high extinction. Star 6, in particular, is surrounded 
by a more substantial region of warm dust and photo-excited material, consistent with lo- 
calized extinction. Star 3 (HD195229), a known B0.2III, is much less red, consistent with 
~1 magnitude of visual extinction, and cons istent with the known variability of extinction 



across Cyg OB2 (IMassey fc Thompson! 1 1 99 ll ). 



We performed aperture photometry on the central star of each bowshock at the mid-IR 
IRAC bandpasses using the SSC basic calibrated data frames, a 10-pixel circular aperture 
for which the aperture correction is minimal, and a much larger annulus to measure the 
diffuse background levels. Most of the stars were saturated at IRAC [3.6] and [4.5] bands in 
the 10.4 s exposures, so the high-dynamic-range (HDR) 0.4 s exposures were used. Table [1] 
lists these fluxes and their uncertainties. 

Figure [TBI shows the spectral energy distributions of the bowshock candidates and central 
stars at the 2MASS, IRAC and two MIPS bandpasses. Asterisks denote photometry of the 
bowshock alone, while diamonds denote the star plus bowshock. The bowshocks are not 
detected at 2MASS JHK bands and are somtimes undetected in the IRAC bandpasses. In 
only three cases (stars 6, 8, and 9) are pointlike sources detected at [24], indicative of dust 



3 http:/ /ssc. spitzer.caltech.edu/irac/calib/ ; 



http: / /ssc. spitzer.caltech.edu/ mips 
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in the immediate vicinity of the stars. Uncertainties are always smaller than the plotted 
points. The solid curves are blackbody fits to the JHK photometry, while the dotted curves 
are fits to the IRAC and MIPS data using scaled iDraine fc Lil (120071 ) dust emission models. 
In Sections 3.4 and 4 we discuss the details of the fitting, the models, and model results for 
each object individually. 



3.3. Bowshock Candidate Morphologies 

The shape of a momentum-driven bowshock is a universal function that scales as the 
"st andoff distance" , -Bp , between the s t ar and the bows h ock a p sis, as s hown analytically 
bv lBaranov et all Jl971): bysonl fjl975h : Ivan Buren et all Jl99p[ l: IwilkirJ Jl996h . Using the 
formulation of IWilkinl (119961 ) Equation 9, we generated synthetic bowshock "images", ap- 



plying arbitrary size, rotation, and intensity scale factors to facilitate comparison with the 
data. Figure [T41 shows a selected subset of bowshock candidates (1, 2, 3, 5, 6, and 7) at [8.0] 
(green), [24] (red) and the scaled synthetic bowshock (blue). The simulated bowshocks are 
all shown at zero degrees inclination, while the inclinations of the bowshock candidates are 
unknown. 

The morphologies of candidates 5 and 7 (central stars 09V and 05V) show striking 
similarity to the theoretical shape. The agreement is best at [24] where the bowshock can 
be traced to larger radial distances than at [8.0] where the bowshock is less extended and 
exhibits a more irregular morphology. For object 2, there is also good agreement, although 
the bowshock itself is faint and difficult to see because of the large dynamic range in the 
images. For objects 1, 3, and 6 the morphological agreement is considerably less good. We 
take this as evidence for a physical origin other than a bowshock. The morphologies of the 
other objects 4, 8, 9, and 10 (not pictured) also differ from the canonical bowshock shape, 
appearing more irregular and clumpy. We interpret the good agreement between the model 
shape and objects 2, 5, and 7 as support for the momentum-driven bowshock nature of these 
objects. 



3.4. Individual Objects 

3.4.1. Bowshock Candidate 1 and Star 1 (GSC03161-01188) 
Bowshock candidate 1, pictured in Figured is associated with the V=13.3 star GSC03161 - 



01188 from the Guide Star Catalog, also known as #1288 from the work of Reddish et al. 
Jl966h . Figure [15] (bottom) shows the WIRO spectrum of GSC03161-01188. We classify 
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G SC03161-01188 as 9V on the basis of the Hel/Hell line ratios and by comparison to 



Jacoby et al.l (119841 ) spectral library 09V star (top). We compute a spectrophotometric 



distance for the K=8.21 star GSC03161-01188 by adopting a visual extinction of 4.5 mag- 
nitudes from Figure d2] (Ak = 0.51 mag), an intrinsic col or V-K=— 0.85, and assuming an 
absolute magnitude for an 09V star of My = 4.05 ± 0.2 (IMartins et al.ll2005l ). This yields 
a spectrophotometric distance of 1500 ± 150 pc for GSC03161-01188, consistent with the 
adop ted distance to Cyg OB2, ~ 1.6 kpc (jHansonll2003l ; iMassey fc Thompson! 1 199 ll ; ISchulte 
1958t h 



Cross correlati ng the spectrum o : 



GSC 03161-01188 with a synth etic model atm o spher e 



09V spectrum from lLanz fc Hubenyl (120031 ). as described more fully in lKiminki et al.l (120071 ). 
yields a heliocentric radial velocity of —17 ± 6 km s _1 . This velocity is consistent with, but 
more negati ye than, the a v erage radial velocities of Cyg OB2 members, —10 ± 6 km s _1 
surveyed by lKiminki et al.l (120071 ) . 



The spectral energy distribution in Figure [13] shows that this nebula is detected only at 
[24] and [70]. Stellar flux dominates the 2MASS and IRAC bandpasses. The solid line shows 
the SEP of an Q9.5V star, approximated as a blackbody with T e jj=30,000 K and R=8 R Q , 
( IMartins et al.ll2005l ) which, at a distance of 1.6 kpc with Ay = 4.5 mag, provides a good fit 
to the near-IR data. This is consistent with our spectral classification based on Figure [T5l 

Given the lack of similarity to the theoretical shape and the modest radial velocity of the 
central star, this nebula is unlikely to be a bowshock. We assign this object an ambiguous 
classification. Possibilities include a dusty asymmetric bubble or a partial shell. 



3.4.2. Bowshock Candidate 2 and Star 2 (G80. 9020+0. 9828) 



Figure [2] shows bowshock candidate 2 and Star 2 (G80.9020+0.9828). Figure M (bot- 
tom) shows the WIRO spectrum of this star a nd a comparison B2 dwarf equivalent from the 
model atmospheres of lLanz &: Hubenyl (120031 ) (top). The relative strengths of the H I 6563 
A and the He I 6678 A and He I 5876 A lines lead us to classify this star as B2V, plus or 
minus several spectral subtypes, owing to the low signal-to-noise spectrum. Its heliocentric 
radial velocity is — 12 ± 15 km s -1 , consis tent with the range of radial velocities observed for 
Cyg OB2 members iKiminki et al.l (120071 ). We compute a spectroph otometric distance fr om 
the 2MASS photometry (K=10.5) by adopting a visual extinction (ICardelli et al.lll989l ) of 
4.5 magnitudes (Ak = 0.51), an absolute magnitude for an B2V star of My = —2.4 ±0.4, an 
intrinsic V-K color of -0.4. This yields a spectrophotometric distance of 2500 ± 500 pc, where 
the error is driven primarily by the uncertainty on the spectral type/absolute magnitude. 
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A high-surface-brightness apsis lies 9" (0.07 pc projected separation) from the star toward 
the upper left in Figure [2J and there is also a more extended region of 8 micron emission 
preceding the apsis to the upper left. The [24] lima-bean shaped nebulosity and the arc-like 
structure seen at [8.0] do not share a common symmetry axis, so it is difficult to ascertain 
which wavelength best traces the putative bowshock. 

Figure [13] shows that Bowshock 2 is a strong detection at the IRAC and MIPS band- 
passes. The relative faintness of the bowshock in [4.5] relative to the other three IRAC 
bandpasses is consistent with strong PAH emission in the shock. The central star dominates 
the 2MASS bandpasses. Star 2 is fainter than Star 1, and is well fit by a B2V-B3V star 
(T e //=20,000 K, R=4 R Q ) at a distance of 1.6 kpc with an extinction of Ay = 4.5 mag (solid 
line). 

Although this candidate bowshock nebula is faint and difficult to represent owing to a 
large dynamic range in the images, the morphological agreement with the theoretical shape 
is reasonably good. The apsis of this nebula is bright at [8.0] and [24] and may indicate 
an interaction or illumination of a of a nearby cloud. We retain this object as a probable 
showshock. 



3.4.3. Bowshock Candidate 3 and Star 3 (HD195229) 

Bowshock candidate 3 appears to be generated by its attendant star, HD195229, listed 
as B0.2III in SIMBAD. Figures [T7| and [18] show our WIRO and WIYN spectra of this star, 
from which we confirm a B0.2III spectral type. From these spectra we find a heliocentric 
radial velocity of — 3±2 km s -1 , slightly less negative than other massive stars near Cyg OB2. 
The luminosity for a BO. 2 giant is somewhat uncertain because the evolution ary tracks for 



mass ive giant and supergiant stars are nearly vertical in an H-R diagram (jMarigo et al. 



2008f ) . Assuming My = -5 ± 0.5, V-K=-0.8 and Ay = 1.0 from Figured for this K = 7.3 
star leads to a spectrophotometric distance of 1800 ± 500 pc, consistent with the distances 
of massive stars in Cyg OB2. 

HD195229 is also known as a hig h- proper- motion star w ith ura = 0.14 ±0.55 mas yr -1 , 



and [iDec — —2.89 ± 0.47 mas yr 1 (jPerryman et al.lll997l ). Its Hipparcos proper motion 



corresponds to a velocity of 28±8 km s _1 in the plane of the sky at 2000 pc. Combining 
this with the star's radial velocity gives a space velocity of 32±8 km s _1 , placing it above 
the nominal 30 km s _1 threshold for runaway status. The yellow bar in Figure [3] shows the 
implied sky motion (from left to right) over 10,000 years and has position angle uncertainty 
of about 15 degrees. The orientation of the candidate bowshock implies a velocity vector 
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toward the lower left of Figure [3] while the proper- motion data imply a velocity that is toward 
the lower right. 

Figure [T3] shows that Bowshock 3 is detected with certainty only at [24] and longer 
wavelengths, with a marginal detection at [3.6] and [4.5]. The stellar SED is well fit by an 
early B giant (T e //=25,000 K, R=10 R®), consistent with its spectral type, at a distance of 
1.6 kpc with an extinction Ay = 1.0 mag. Curiously, HD 195229 has a much lower extinction 
than the rest of the stars in Cyg OB2 (4-6 mag), requiring either that this star is seen 
through a local minimum in the obscuring dust, or that HD 195229 is actually on the near 
side of a heavy veil of extinction that enshrouds most of Cygnus X. 

The lack of similarity between the [24] appearance and the theoretical shape in Fig- 
ure d3l coupled with the misalignment between the proper motion and putative bowshock 
morphology, casts doubt on a bowshock interpretation. We categorize the nature of this 
nebula as ambiguous. 



3.44. Bowshock Candidate 4 and Star 4 (Gl 9. 4111+1.21 03) 



Bowshock candidate 4 resembles Bowshock 1, having no detectable emission at the 
IRAC bandpasses and a lima-bean morphology at [24]. Note the edge-brightened cloud to 
the right of the bowshock, suggesting illumination by the (presumed) early-B star powering 
the nebula. The central star is the faintest in our sample at optical and infrared wavelengths 
and Figure [121 suggests a very large extinction of Ay > 12. Its K-band magnitude is similar 
to the central star of Bowshock 2, so it is probable that this is also an early B star, if 
located at a similar distance. There is no literature identification for this star, so we give it 
a designation according to its Galactic coordinate, G79. 4171+1. 2703. 

We note that the putative bowshock nebula points away from Cyg OB2 and in the direc- 
tion of the H II region IRAS20264+4042, located toward the upper righ t from Star 4 in Fig- 
ure [D This star forming region is also identified as the radio source DR7 (IDownes &: Rinehart 
1966h and the host of the ultra-compact H II region G79.320+1.313 fcurtz et alill994f ). 
While it is possible that outflows from this SF region are responsible for the high relative ve- 
locities leading t o the produc t ion of the putative bowshock (as shown for bowshocks in M17 
and RCW 49 by lPovich et al.l (120081 )). a careful search for other similarly oriented bowshocks 
around IRAS20264+4042 revealed none. 



Odenwald fc Schwartz! fll993J ) report a 12 CO detection for IRAS20264+4042 a t V hf ,R 



—42 k m s l , corresponding to a kinema tic distance of 6-7 k pc i nferred from the 



(119851 ) Galactic rotation curve. However, iDutra &: Bical (1200 ll ) and lLe Duigou fc Knodlseder 



Clemens 



13 



( 120021 ) suggest distances of 1.1 and 1.6 kpc respectively, based on near infrared photometry. 
At such a distance, the K=11.3 magnitude of the central star, coupled with an implied 
extinction of Ax ~1.5 mag would correspond to an early- to mid-B star. At the the larger 
(but, we contend, less probable) distance of 6.5 kpc, the apparent magnitude would be 
consistent with an 09V or similar. Given the abundance of early-type stars at the nearer 
distance of Cyg OB2, we prefer the near distance and the corresponding early-to-mid B 
spectral type. 

Figure [TBI shows that Bowshock Candidate 4 is detected only at [24] and [70]. Although 
the stellar spectral type is not known with certainty, the stellar SED is well fit by an ~B3 
dwarf (T e ^j=20,000 K, R=4 R@) at a distance of 1.6 kpc with a very high extinction of 
Ay = 16 mag. The stellar photometry exhibits an excess at [8.0], suggesting the presence 
of hot dust or possibly PAHs associated with the star. Lack of similarity to the theoretical 
bowshock shape leads us to classify this object as ambiguous. 



3.4-5. Bowshock Candidate 5 and Star 5 (A10) 
Figure [5] shows Bowshock Candidate 5 to be a stunning combination of emission at 



the I RAC and MIPS bandpasses. The central star is A10 in the notation of IComeron et al. 



( 120021 ). Our new spectrum of A10, shown in Figure [IS is consistent with a late O star, 
although the uncertainty is several spectral subclasses. Its heliocentric radial velocity is 
10±10 km s _1 , consistent with other Cyg OB2 members. Its JHK colors and magnitudes 
are similar to the central star of Bowshock Candidate 1, suggesting a similar reddening and 
late O spectral type if at the distance of Cyg OB2. The angular distance between the star 
and bowshock apsis is ~25" or 0.19 pc projected separation. 

Figure d shows that the 2 MASS data are well fit by an 09V (T e// =31,000 K, R=8 
Rq) at a distance of 1.6 kpc with an extinction of Ay = 7.5 mag. The nebula dominates the 
SED at the IRAC and MIPS bandpasses. 

Figure HH shows that the [8.0] and [24] appearance closely resembles the theoretical 
bowshock shape. The orientation is consistent with either motion toward Cyg OB2 or in- 
teraction with a flow of material emanating from Cyg OB2. We classify this object as a 
probable bowshock. 
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3.4-6. Bowshock Candidate 6 and Star 6 (G80. 7657+0. 4966) 

Bowshock candidate 6 exhibits strong emission at both IRAC and MIPS bandpasses 
in Figure [6j It lies within a complex of diffuse emission that may constitute a small H II 
region. The central star is extremely red, and Figure [12] suggests an extinction approaching 
Ay = 18. Accordingly, there is no literature identification or spectral type available, so we 
designate it by its Galactic coordinates, G80. 7657+0.4963. 

Though the nebular component of object 6 is strong at the IRAC and MIPS bandpasses, 
the central star is the faintest in the sample at IRAC bandpasses. While there is no spectral 
information on Star 6, if we assume it is at 1.6 kpc the 2MASS SED is well fit by a mid-B 
dwarf (Te//=15,000 K, R=3 R & ) with a very high extinction of Ay =19 mag. 

Star 6 is one of only three detected at [24] , suggesting the presence of a dusty shell or 
disk associated with the star, since we would not expect to detect the stellar photospheres 
at [24] or [70] . Figure [12] shows that the flux from the star at [24] is greater than the nebular 
flux and continues to rise beyond 24 microns. These characteristics suggest that Star 6 is a 
pre-main-sequence object illuminating a rim of material that points toward Cyg OB2, the 
likely source of a high-velocity wind. 

Figure [7] shows the [8.0] image (green), [24] image (red), and 12 CO (1-0) molecular 
line ma p (blue and c ontours) integrated between LSR velocities -7 - -3 km s _1 and 9-14 



km s~ 1 (IBrunt II2009I ). This figure shows that Bowshock 6 is located along a tendril of CO 
(at velocities -7 - -3 km s _1 ) pointing toward Cyg OB2. It lies adjacent to a prominant 
gasseous "pillar" seen at [8.0] and in CO. Th e star forming region I RAS 20343+4129, which 



contains several pre-main-sequence B stars (ICampbell et al.ll2008l ). constitues the head of 



this pillar and is coincident with a concentration of CO at velocities 9-14 km s _1 . This 
pillar also points toward Cyg OB2, consistent with it being eroded by radiant and mechanical 
luminosity of the massive stars therein. The similarity between the contents (young early B 
stars) and morphology (pillar/CO filament) of the Star 6 region and IRAS 20343+4129 raises 
the intriguing possibility that these two represent different phases, and perhaps different 
viewing angles of the same phenomenon: massive star formation triggered in the heads of 
gas columns being eroded from the outside. The lack of a strong CO concentration and 
visible pillar of photo-dissociated material near Star 6, coupled with the presence of the 
bowshock pointing toward Cyg OB2 suggests that the region containing Star 6 may be at 
a more evolved stage compared to IRAS 20343+4129. It may also be that the stars within 
IRAS 20343+4129 are more massive and, thereby, generate a more luminous H II region. 
We classify the central star of Bowshock Candidate 6 as a probable young stellar object. 
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3.4-7. Bowshock Candidate 7 and Star 7 (A37) 



Bowshock candidate 7 has the largest angular extent of any of our targets, covering 
nearly 5' in Figure [HJ It i s detected only at [24 ] and [70]. The central star appears to 
be A37 in the notation of Comeron et al. iHanson] (l2003h determines a spectral 

type of 05V((f)), making it the most energetic and ma ssive in our sample. Th e prese nce 
of a bowshock associated with the this star is noted by iGvaramadze fc Bomand (120081 ) on 
the basis of images from the Midcourse Space Experiment mission. The angular distance 
between the star and bowshock apsis is the largest in our sample, ~70", or 0.53 pc projected 
separation. 



The SED of Star 7 in Figure E3] is consistent with an 05V star (T e// =40,000 K, R=ll 
R & ) and Ay = 5 mag if placed at a distance of 2.1 kpc, slightly larger than the nominal 
distance to Cyg OB2, al though the line-of-sight dept h along the Cygnus-X complex is likely 
to be ~few hundred pc ( IQdenwald fc Schwartzlll993l ). 



The excellent agreement between the [24] morphology and the theoretical shape leads us 
to classify this as a probable bowshock. The orientation of the bowshock indicates a motion 
toward the left in Figure [BJ roughly tangent to the vector toward Cyg OB2. This suggests 
that the origin of A37 lies somewhere outside of Cyg OB2. 



3.4.8. Bowshock Candidate 8 and Star 8 (G77. 5168+1. 9047) 

Bowshock candidate 8 is visible at both IRAC and [70] bandpasses. Figure |9] shows 
that this object lacks the classical arc morphology and symmetry of the other bowshocks, 
warranting an ambiguous classification. It resembles a partial shell or bubble. The WIRO 
spectrum of this star obtained 2009 September 16 is consistent with a B2V - B3V having a 
heliocentric radial velocity of 2±4 km s _1 , similar to other Cyg OB2 members. The apparent 
central star, designated G77. 5168+1. 9047, is relatively bright at J=10.5, and it has 2MASS 
magnitudes and reddening consistent with other bowshock stars and early-type Cyg OB2 
members. 

Figure EH shows that the near-IR SED of Star 8 is well fit by a B3V (T e// =20,000 K, R=4 
Rq) with Ay = 4.5 at a distance of 1.6 kpc. The IRAC and MIPS bands are dominated 
by extra-photospheric emission, consistent with hot dust near the star. Star 8 is one of 
only three stars (Star 6 and Star 9) detected at [24], suggesting the presence of dust in a 
circumstellar disk or envelope. The morphology of the nebula is does not bear a strong 
resemblance to the theoretical bowshock shape. As such, we designate Star 8 as a possible 
young stellar object (YSO). 
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3.4-9. Bowshock Candidate 9 and Star 9 (G76. 8437+0. 1231) 



Bowshock candidate 9 is a compact nebula with detected emission only at [5.8] and 
longward. The central star is clearly seen as a point source in Figure [TTJJ The stellar SED 
of Star 9 is difficult to ascertain, since there is no 2MASS counterpart to the point-like 
source detected at the mid-IR IRAC bands. The mid-IR flux rise steeply with wavelength, 
consistent with a heavily enshrouded source. The [4.5] flux lies below the expectations of a 
simple interpolation between the [3.6] and [5.8] bands. This has beco me a classic s i gnatu re 
of material dominated by emission from large molecules (e.g, PAHs; iDraine fc Lil (120071 )). 
which are largely absent in the [4.5] band. The nebular morphology does not bear a strong 
resemblance to the classic al bowshock shape. These characteristics are consistent with a 
young (class I?) protostar (IKenyon et al.lll993l ). and we designate this object as a probable 
YSO. 



3.4.10. Bowshock 10 and Star 10 (GG77. 0511-0. 6092) 

Bowshock 10, pictured in Figure [TT] is a small, narrow filament with strong emission in 
the IRAC bands, and a low ratio of far-IR to mid-IR flux in the final panel of Figure [TBI 
The stellar source is identified only by its Galactic coordinates, G77.0511-0.6092. Because it 
has similar 2MASS colors and magnitudes as many of the other bowshock stars, we consider 
it to be a probable late O or early B star at the distance of Cyg OB2. 

The near-IR SED of Star 10 is consistent with an early B star (T e jj=25,000 K, R=5 
Rq) with Ay = 6 at a distance of 1.6 kpc. The IRAC points lie in excess of the photosphere, 
consistent with hot dust and emission from PAHs. The morphology of the [8.0] nebula 
appears irregular, while the [24] nebula is similar to the classical bowshock morphology. We 
classify this object as ambigous. 



3.4.II. BD+43°3654 



Comeron fc Pasqualil (120071 ) present proper motions and MSX i nfrared images of the 



Q4If s tar BD+43°3654, originally identified as a probable runaway by Ivan Buren fc McCray 
(119881 ) on the basis of a probable bowshock seen with IRAS0 They compute a distance of 



4 BD+43°3654 lies outside the area of the Spitzer Legacy Survey of the Cygnus-X Complex, so we do not 
discuss this object in detail or present figures. 
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1450 pc and p ropose that BD+43°365 4 was ejected from Cyg OB2 1.6 Myr ago, not long 
after its birth. iGvaramadze &: Bomand (120081 ) propose that BD+43°3654 is a blue straggler, 
formed from the merger of two stars and ejected from Cyg OB2 during a binary-binary 
encounter. They invoke the high proper motion pulsars B2020+28 and B2021+51 as the 
descendants of the other two stellar participants in this 4-body interaction. 

We cross-correlated the WIYN spectrum of BD +43°3654 with an 051 model atmo- 
sphere, which yielded a heliocentric radial velocity of Vlsr=— 66.2±9.4 km s _1 . This is con- 
siderably more negat ive than the Vlsr—— 10 km s _1 mean velocity of O stars in Cyg OB2 
( IKiminki et al.1 120071 ). and therefore BD +43°3654 would qualify as a runaway on the ba- 
sis of its radial velocity al one. Its motion in the plane of the sky is 39.8 ± 9.8 km s _1 
(jComeron fc Pasqualil 120071 ) at the assumed distance. Combining this with our radial veloc- 
ity measurement, we calculate the heliocentric space velocity of BD+43°3654 to be 77 ± 10 
km s _1 . 



3.5. Summary of Classifications and Statistics of Runaways in Cygnus X 



The bowshock candidates presented here comprise a mixed bag of phenomena. Objects 
2, 5, and 7 (and BD +43°3654) are most likely to be genuine bowshocks powered by runaway 
stars on the basis of having early spectral types and nebular morphologies most similar to 
the classical bowshock appearance. However, none of stars 2, 5 or 7 have remarkable radial 
velocities. Objects 1, 3 (the high proper-motion star HD195299), 4, and 10 are ambiguous, 
having nebular morphologies less similar to the classical bowshock shape, but, nevertheless, 
appear to be early-type stars. Objects 6, 8, and 9 have central stars with excess IR emission, 
suggesting circumstellar material. Stars 6 and 9 also have high extinction and/or SEDs that 
rise toward longer wavelengths, consistent their being young stellar objects hosting disks or 
circumstellar envelopes. The nebula associated with Star 6 shows evidence for being the 
head of a gaseous pillar being ablated by radiation from the direction of Cyg OB2. Hence, 
objects 6, 8, and 9 are the least likely to be high velocity runaway stars. 



In a magnitude-selected sample of 195 bright O stars, iGies fc Boltonl (119861 ) found that 
16 (8%) could be classified as runaways either on the basis of radial velocities (6; 3%), proper 
motions (5; 2.5%), or extreme distance from the Galactic Plane ( 5; 2.5 %). By contrast 
the radial velocity survey of 146 Cyg OB2 stars by IKiminki et al.1 (120071 ) did not find any 
runaways. The SST Cygnus X Legacy Survey covers a muc h larger region (~ 2 3 square 
degrees) than the sample of core Cyg OB2 members studied bv lKiminki et al.1 (120071) f~ 0.3 5 



square degrees), which was selected from the optical imaging survey of iMassey &: Thompson 
( )199ll ). The ratio of areas in these two surveys is roughly 65, so one might conclude that the 
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number of early type stars in Cygnus X is as many as 146 x65 = 9500. However, the earliest 
and most massive sta r s app ear to be contained within the core Cyg OB2 region covered by 
Massey &: Thompson! (Il99ll ). so that this is a gross overestimate. We very conservatively 
adopt that the SST Cygnus X Legacy Survey survey inc ludes ~600 O B stars earlier than B2, 



or about four times the number of OB stars surveyed by lKiminki et al.1 (120071 ). In this region 
we find only three candidates for runaways — four if BD +43°3654 is included. At face value, 
this is 3/600 = 0.5%. Of course, high space velocities are a n e cessar y but not a sufficient 
condition for for the formation of bowshocks. Ivan Buren et al.1 (119951 ) found that 58 out of 
188 high velocity stars evince bowshocks at IRAS sensitivity levels, and of these, 25 have well 
resolved structure. We expect that SST would detect even fainter and smaller bowshocks. If 
we assume that 1/4 of runaway stars produce bowshocks visible at SST sensitivity levels and 
that 8% of all OB stars are runaways, we would expect to find on the order of 12 bowshocks 
in the Cygnus X region. Why do we detect so few? 

We identify only three probable bowshock runaways, and perhaps only one of these 
(Bowshock 2) has an orientation suggesting an origin near Cyg OB2. It is certainly possible 
that most of the runaways, like BD +43°3654, could fall outside the survey area if they 
have high space velocities and wereprobable ejected from the Association very early in its 
formation. Stars with a tangential velocity of 30 km s~ x would travel beyond the boundaries 
of the SST Cygnus X Legacy Survey (~26 pc diameter) in about 0.5 Myr — a small fraction 
of the 2-4 Myr age of Cyg OB2. Dynamical ejection scenarios wherein N-body interactions 
in close multiple systems expel stars would predict this — that most of the runaways are 
produced early in the formation of a young massive star cluster and have already traveled 
beyond the immediate cluster vicinity. Supernova ejection scenarios, on the other hand, 
predict that the start of the ejection process is delayed several Myr (until the explosion 
of the most massive stars), and then it continues at a constant or increasing rate as stars 
from lower-mass OB binaries complete their evolution. Although the arguments are mostly 
qualitative, our findings would seem to be more consistent with the dynamical ejection 
scenario. 



4. Bowshocks and Laboratories 



Baranov et al.l (ll97ll ) first articulated the theoretical basis for momentum-driven bow- 



ar wind bow- 



shocks in the context of the solar wind termina tion shock. The theory of ste l 
shocks has since been f urther refined in papers by lvan Buren fc McCravl (1988) .Ivan Buren et al 
dl990h , Iwilkinl Jl996h , IwilkirJ (boooh and simulated numerically bv lComeron fc Kaper! Jl998h. 



A bowshock forms at a radial distance Ro where the momentum flux from the stellar wind 
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material of density p w and velocity V w is equal to the momentum flux from the ambient 
medium of density p a moving at V a relative to the star. The relative velocity may be the 
result either of a "runaway" star moving at high speed or a bulk flow of the ISM, such as an 
outflow from a young star-forming region. Hence, 



1 2 1 2 

7>PwV w = ~p a V a . 

The density of the stellar wind can be expressed as 



_ M w 

Pw ~ 4irRlV w ' {Z) 

where M w is the stellar mass loss rate. It follows that the standoff distance of the bowshock 
is 



R 



V W M W 



(3) 



Rq can be measured directly from the bowshock images, assuming that all of the stars 
are located at the adopted distance of Cyg OB2. Technically, only the projected separation, 
Rq cos i, can be measured, but in order to observe a bowshock morphology, the viewing 
angle cannot be far from % = 0° so that cos i is not far from unity. Adopting a mean ISM 
gas mass per H atom of p = 2.3 x 10~ 24 g means that p a can be expressed in terms of the 
ambient number density: p a = pn a . Stellar wind speeds and mass loss rates h ave proved 



difficult to measure and m ay be even more uncertain than previously believed (IPuls et al. 



20081 ; iFullerton et al.l 120061 ). As an estimate of the ma ss loss rates and win d velocities for 
stars in our sample, we adopt the values tabulated by iMokiem et al.l (120071 ) for early type 
stars. Equation 3 can then be rearranged to yield the relative motion btween the star and 
ambient medium, 



v - - \S-« 1/2 • <4) 

In practice, it is difficult to calculate the star's velocity, V a , relative to the surrounding 

medium because the quantities on the right hand side are poorly known, especially the 

density of the ambient medium, n a . By adopting published estimates for the stellar mass 

i ii 

loss rates and wind velocities, we can at least make an estimate of the product V a n a , 
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similar to the approach of IPovich et al.l (120081) . In kee ping with the precedent established in 
van Buren fc McCray ( 1988 ) and Povich et al. J2OO8 ). we write the mass loss rate in units 
of 10 -6 M Q yr _1 as M w _ 6 , the stellar wind velocity in units of 10 8 cm s _1 as Kg and the 



ambient medium density in units of 10 3 cm 3 as n a 3. The relative star-ISM velocity is then, 



V a = 1.5(— )'\V w ,s M^-e 
pc 



l/2 n -l/2 



a,3 



[km s~ 



(5) 



Given the uncertainties in stellar mass loss rates and wind velocities for stars later than 
about BO, we compute and tabulate derived parameters only for the five objects in our 
sample with the most reliable spectroscopic types, objects 1, 3, 5, 7, and 11 (BD+43°3654). 
Table [3] tabula t es the adopted stellar wind velocities, mas s loss r ates (from the works o f 
Repolust et aD (hooj): iMokiem et al.l J2005h: iMartins et all (bo05[ ): ICrowther et all (bood ). 
as summarized in lMokiem et al.l (120071 )). standoff distances, Rq cos i, along with the product 



1/2 

V a n a ' 3 . The standoff distances of 0.11 - 1.8 cos i pc are consistent with those found for 



1/2 

the bowshocks in IPovich et al.l (120081 ). The derived values for V„ n„ q lie in the ran ge 0.2 
to as much as 30 — broadly consistent with the values from the (IPovich et al.l 120081 ) stars. 
For typical ISM densities of n aj 3=0.1 - 1.0, the implied stellar speeds are fewxlO km s _1 , 
consistent with the expected speeds of runaway stars. 



4.1. A Novel Measure of Stellar Mass Loss Rates 



Because stellar mass loss rates are pr obably the m o st unc ertain parameter in Equation 
5 above, we invert the approach taken by IPovich et al.l (120081 ) and exploit the phenomenon 
to yield an estimate of the mass loss rates for the central stars. Inverting Equation 5 yields 
the mass loss rate of the star as a function of the relative motion between the str and ISM, 
the ambient density, the stellar wind speed, and the standoff distance, 



M„ 



0.67[i? (pc)] 2 [K(Ws)]V3 



V, 



w,8 



(6) 



Rq is readily observed from the bowshock images, assuming a 1.6 kpc distance to Cygnus X 
and assumi ng cosi g 1. We assu me a minumum star-ISM speed of V a =30 k m s~ l for run 
away stars (jGies fc Boltonl Il986l ). and we adopt the stellar wind speeds of IMokiem et al. 
( 120071 ) which, for our sample of 05 - B0 stars, range from 3000 km s _1 to ~800 km s _1 
(V Wt 8=3 to <0.8). The ambient number density, n a>3 , is difficult to measure directly, but 
may be estimated from n s 3 , the post-shock density wi thin the luminous bowsh ock nebu- 
lae. For strong, highly supersonic shocks, n s> 3 ~ 4n ai 3 (ILandau fc Lifshitzlll987l ). Without 
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knowing the Mach number of each individual case, we adopt n s ^ = 2n a ^, and consider the 
uncertainties to be factors of 2-4. We estimate n Sj3 by fi tting the bowshock spectral energy 
distributions in Figure [12] with the dust/PAH models of iDraine &: Lil (120071 ) which give the 
dust emissivity per H nucleon as a function of the radiation field intensity, parameterized as 
a multiple, U, time s the mean interstellar radiation field of £7/5^=2.3 xlO -2 erg s -1 cm -2 



( jMathis et al.lll983l ). We make the motivated assumption that the bowshock luminosity 
stems from reprocessed photo ns from the cen tral star rather than thermalized mechanical 
energy from the stellar wind ( Benjamin 20091 ). The appropriate choice of model radiation 
field intensity, U, is then the ratio of the flux from the central star at the standoff distance 
over the mean interstellar field, estimated as, 



U 



Rl 



u ISRF > 



(7) 



where R* is the stellar radius, T* the effective stellar temperature, and a the usual Stefan- 
Botlzmann constant. The postshock number density can then be estimated as 



n 



s,3 



where V is the bowshock volume which we estimate by approximating the bowshocks as 
hollow cones at a known distance, D, and N = L Q b s /Lo, the ratio of the observed lum i nosity 
at a given bandpass to the theoretical luminosity per nucleon in the IDraine fc Lil (120071 ) 
models. Eff ectively, n.,3 is fou nd empirically because it is the scale factor by which we 
multiply the IDraine fc Lil (120071 ) models to match the observed SED. Because V oc D 3 and 
the inferred luminosity L & s oc D 2 , our estimate for n s ^ varies as D~ l , inversely with the 
adopted distance. 



Draine &: Lil (120071 ) parameterize the model spectral energy distributions in terms of a 
maximum and minimum radiation field indicent on the dust, U max and U m i n . While dust in 
the typical ISM is exposed to a range of intensities, we expect that material in the bowshocks 
is irradiated by a single intensity dictated by the stellar luminosity and standoff distance. 
Therefore, we adopt the models with U ma x. = U min = U, where U is chosen according to 
Equation 7. Figure 13 of IDraine fc Lil (120071 ) shows that the flux ratio [24]/ [70] is sensitive 
to U over the range 1 < U < 10 7 , and we make use of this sensitivity in fitting the bowshock 
SEDs below. We adopt the models appropriate to Milky Way dust and leave the PAH 
fraction, qpAH, as a free parameter. Larger values of (Jpah raise the flux in the IRAC 
bandpasses where PAH emission is strong, but have no effect on the SED at [24] and [70]. 



Dotted lines in Figure [13] show illustrative dust models that provide reasonable fits to 
the mid- and far-IR data. The slope of the SED between [24] and [70] allows us to guess 
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the appropriate value for U without any prior knoweldge of the the stellar luminosity or Rq. 
Interestingly, the values for U in each case calculated from Equation 7 are generally within 
factors of 2-3 of the initial guess. This correspondance serves as a kind of corroborating 
evidence that the bowshocks are, in fact, powered by the radiant energy from their central 
stars. The fits of the models to the data in Figure [13] are done by eye to maximize the 
agreement at [24] and [70] by adjusting n Sj3 , and then the parameter qpAH is varied to 
achieve a best fit at the IRAC bands. The models provide reasonably good fits to the data 
for objects 3, 5, & 8. In other cases, the overall fit is poor at the IRAC bandpasses because 
the model PAH features are too weak, even for the maximum model qpAH, or too weak, even 
for the minimum model Qpah- We take these discrepancies as indications that PAHs may 
be destroyed in the bowshocks (former case; objects 1, 4, & 7), and that PAH excitation 
may be elevated and dominate the SED in some objects (latter case; objects 2 & 6). The 
adopted values for U and n s ^ appear in each panel. 

Armed with estimates of the ambient number densities, n a ^ = l/2n Sj 3, we proceed to 
calculate mass loss rates, M w -q, using the method outlined above. Again, we consider only 
the four objects with the most reliable spectroscopic types where the bowshock nature is 
more likely to be correct. Table H] summarizes the relevant parameters adopted for each of 
objects (2, 5, 7, 11). Derived mass loss rates vary from < 0.18 x 10~ 6 M Q yr _1 for the early 
B star (Star 2) to ~ 1 x 10~ 6 M yr" 1 for the 09V star (Star 5), to ~ 3 x 10~ 6 M Q yr" 1 for 
the 05V star (A37). The 04If star BD+43°3654 stands out, having an implied mass loss 
rate of 160 xlO -6 M Q yr -1 . This large value results from the rather large standoff dis tance, 



Rq = 1.4 pc, implied by the MSX infrared image from IComeron Sz Pasqualil (120071 ). We 
note that we have no direct measure of the ambient ISM density for this object, so we 
have assumed the average value from the other four early type bowshock stars. Results 
for BD+43°3654 should be regarded as particularly uncertain. These results are broadly 



consis t ent with the mean M v alues for mid- to late-0 stars estimated by lLamer fc Leitherer 



compilation of iMokiem et al.l (120071 ). 



( 119931 ); iFullerton et al.l (120061 ) but a factor of several or more larger than the more recent 



We consider the uncertainties on mass loss rates to be factors of 2-4, dominated by the 
uncertainties on n a S , which are of the same magnitude. The other parameters that appear 
in Equation 6 are likely known to ~30%, so that uncertainties in the standoff distance, which 
enter as R^ and V^, do not have a major impact on the result. We present the calculation 
of mass loss rates in this section as an outline of a novel method without claiming that this 
approach is in any sense superior to other well-established techniques. 
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5. Conclusion 

We have identified ten candidate bowshock nebulae and their central stars in the region 
surrounding Cyg OB2 on the basis of mid-IR morphologies. We provide the first spectral 
types and radial velocities for several stars (objects 1, 2, 5, & 8) from this survey, showing 
that they are late-0 or early B dwarfs. We find HD 195229 to be a likely runaway on the 
basis of its space velocity calculated from both radial velocity measurements and Hipparcos 
proper motions, but the velocity vectors implied by the proper motion data and bowshock 
orientation are nearly orthogonal. We measure the radial velocity of the suspected runaway 
04If star BD+43°3654 to be Vheiio = — 66±4 km s" 1 , supporting its runaway status. Based on 
the morphologies, spectral types, and spectral energy distributions, we identify three objects 
(6, 8, 9) as probable B-star young stellar objects with circumstellar material. Objects 2, 
5, and 7 are probable bowshocks based on similarity to the theoretical shape. These stars, 
however, have modest radial velocities, consistent with those of other Cyg OB2 members, 
indicating that their motions must be predominantly tangential to the line of sight. Without 
proper motion data, we are not able to confirm a runaway status or suggest an origin for these 
stars. The physical nature of the nebulae around stars 1, 3, 4, and 10 remains ambiguous 
with the current data. These nebulae may be partial shells, partial bubbles, or rims of nearby 
clouds illuminated by the early-type stars. As part of our analysis, we have also proposed 
a novel method for determining the mass loss rates for runaway OB stars using bowshock 
physics that yield reasonable, if slightly larger values for M compared to more traditional 
methods. 
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Table 1. Stellar Sources Associated with Bowshock Candidates 
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Table 2. Photometry of Bowshock Candidates 
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Note. — Uncertainties are estimated at 20% in all cases except nondetections where la upper limits are given. 
a High and variable background precludes a reliable measurement at [24] for this source. 



Table 3. Calculated Bowshock Parameters for Most Reliable Objects 
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a Based on spectroscopic identification. 

b Based on the images of IComeron fc Pasqualil | |2007|) and adopted 1.6 kpc distance. 



Table 4. Derived Mass Loss Rates for Most Reliable Objects 
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othe ther stars. 



- 30 




Fig. 1. — Spitzer three-color image in Galactic coordinates of the Cygnus-X region centered 
on Cygnus OB2 with [4.5] in blue, [8.0] in green, and [24] in red. Seven of the ten bowshock 
candidates are labeled with their bowshock number and identifications for the central stars 
from Table [U Several other prominent IRAS infrared sources are labeled, some of which are 
the heads of gaseous pillars pointing toward Cyg OB2, such as IRAS20343+4129. 
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Fig. 2. — Spitzer three-color image of Bowshock Candidates 1 (GSC03161-01188) and 2 
(G80. 9020+0. 9828) with the same color scheme as Figure HJ These two objects present 
a striking contrast as Bowshock 1 is seen only at 24 microns and longer wavelengths while 
Bowshock 2 has strong PAH emission in the IRAC bands and is weaker at longer wavelengths. 
The arrow shows the direction toward the core of Cyg OB2. 
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Fig. 3. — Three-color image as in Figure [T] showing the bowshock candidate around Star 
3 (HD195229). The yellow line indicates HD195229's proper Hipparcos motion, indicating 
that there is a misalignment between the velocity vector inferred from proper motion data 
and the vector inferred from the putative bowshock morphology. The uncertainty on the 
position angle of the proper motion is ~15°. 
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Fig. 4 



- Three-color image of Bowshock Candidate 4, as in Figure [TJ 



Fig. 5. — Three-color image of Bowshock Candidate 5, as in Figure [TJ 
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Fig. 6 



- Three-color image of Bowshock Candidate 6, as in Figured! 
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Fig. 7. — Three-color image of Bowshock Candidate 6, with [8.0] in green, [24] in red, and 
the 12 CO (1-0) map between LSR velocities -7 - -4 and 9-14 km s _1 in blue and contours. 
Star 6 and Bowshock Candidate 6 lie adjacent to a gasseous pillar containing several massive 
young protostellar objects at its head (IRAS 20343+4129). 
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Fig. 8 



- Three-color image of Bowshock Candidate 7, as in Figure [TJ 
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Fig. 11. — Three-color image of Bowshock Candidate 10, as in Figure [U 
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Fig. 12. — A JHK color-color diagram showing the sample stars (diamonds) in relation to the 
fiducial main sequence and supergiant sequence (solid line and dashed line). The reddening 
vector shows the equivalent of 5 magnitudes of visual extinction based on the ICardelli et al. 
(119891 ) reddening curve. 
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Fig. 13. — Infrared spectral energy distributions of the sample nebulae based on the pho- 
tometry from Table [2j 
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Fig. 14. — Images of selected bowshocks with [24] in red, [8. 0] in gr e en, a nd a simulated 
bowshock, based on the analytical shape parameterization of IWirkinl (119961 ). in blue. The 
angular scales and rotations are arbitrarily adjusted to facilitate comparison between the 
data and models. 
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Fig. 15. — Standard 09V spectrum from lJacoby et al.l (11984] ) (top) and WIRO spectrum of 
Star 1 (GSC03161-01188) (bottom). The relative strengths of the H I 6563 A, He I 6678 A 
and He II 5411 A lines lead us to classify GSC03161-01188 as a 09V. IS denotes interstellar 
bands. 
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Fig. 16. — A model atmosphere B2V equivalent spectrum from lLanz fc Hubenyl (120031 ) (top) 
smoothed to the resolution of the data and WIRO spectrum of Star 2 (bottom). The relative 
strengths of the H I 6563 Aand the He I 6678 A and He I 5876 A lines lead us to classify 
this star as B2V, plus or minus two spectral subclasses. IS denotes interstellar bands. 
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Fig. 17. — Standard B0.2III spectrum from iJacoby et al.l (1 19841 ) (top) and WIRO spectrum 
of Star 3 (HD195229) (bottom). The WIYN spectrum of Star 3 in Figure [TS] was used to 
determine this star's spectral type. 
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Fig. 18.— Standard B0.2III spectrum of HD 108639 flWalborn fc Fitzpatrickl 119901 ) and 
WIYN spectrum of Star 3 (HD195229) (bottom). The strengths of the Si IV 4089, Si IV 
4116 A and C III 4069 A lines identify the star as B0.2III. 
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Fig. 1 9. — WIRO spectrum of Star 5, also known as A10 in the notation of IComeron et al. 
(120021 ). The spectrum is consistent with a late O star, although we estimate the uncertainty 
to be several spectral subclasses owing to the low signal-to-noise ratio of the spectrum. 
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Fig. 20. — WIRO spectrum of BD+43°3654. Comparison with a standard 05If from 
Jacoby et al.l (11984) fthe earliest O supergian t in this library) agrees with the 04If clas- 
sification given by lComeron fc Pasquali ( 2007 ). 



